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Q ABSTRACT 

v, The relic abundance of a particle species that was once in thermal equilibrium 

3 in the expanding Universe depends upon a competition between the annihilation 
rate of the species and the expansion rate of the Universe. Assuming that the 
Universe is radiation dominated at early times the relic abundance is easy to 
compute and well known. At times earlier than about 1 sec after the bang theie 
is little or no evidence that the Universe had to be radiation dominated, although 
that is the simplest — and standard — assumption. Because early-Uni verse relics 
are of such importance both to particle physics and to cosmology, we consider 
in detail three nonstandard possibilities for the Universe at the time a species 
abundance froze in: energy density dominated by shear (i. e., anisotropic ex- 
pansion), energy density dominated by some other nonrelativistic species, and 
energy density dominated by the kinetic energy of the scalar field that sets the 
gravitational constant in a Brans-Dicke- Jordan cosmological model. In the sec- 
ond case the relic abundance is less than the standard value, while in the other 
two cases it can be enhanced by a significant factor. We also mention two other 
more exotic possibilities for enhancing the relic abundance of a species- a larger 
value of Newton’s constant at early times (e.#., as might occur in superstring 
or Kaluza-Klein theories) or a component of the energy density at early times 
with a very stiff equation of state (p > p/ 3), e.g a scalar field <f> with potential 
V(<j>) = fi\<f)\ n with n > 4. Our results have implications for dark matter searches 
and searches for particle relics in general. 
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I. INTRODUCTION 


The existence of the cosmic microwave background radiation (CMBR) with 
temperature 2.74 K is very strong evidence that the Universe was both radiation- 
dominated [for times earlier than tgQ ~ 4.4 x 10 10 (flo^ 2 ) -2 sec] and very hot 
[T ~ MeV(t/sec) -1 / 2 ] at early times. 1 Because the temperatures reached early 
on were so high there is every reason to believe that essentially all the known 
particle species and possibly other species yet to be discovered were present in 
great numbers. If equilibrium thermodynamics were the entire story these facts 
would be of little interest, as today the equilibrium abundance of any massive 
particle species would be exponentially small, proportional to exp(— m/T). How- 
ever, it has long been realized that due to the expansion of the Universe the 
actual abundance of a stable particle species cannot track its equilibrium value 
forever, and depending upon the strength of its interactions, the abundance per 
comoving volume eventually ceases to decrease and freezes in at some constant 
value. “Freeze in” of the particle’s abundance occurs when the annihilation rate 
can no longer keep pace with the expansion rate of the Universe: Roughly, the 
abundance ceases to decrease when the annihilation rate falls below the exp an - 
sion rate — when annihilations are said to freeze out. (The reactions that regulate 
the number of a particle species are pair production and annihilation; the pair 
creation rate is related to the annihilation rate by detailed balance, or time- 
reversal invariance.) Moreover, the weak shall dominate — the relic abundance of 
a particle species is inversely proportional to its annihilation cross section. 

Calculating the relic abundance of a particle species that was once in thermal 
equilibrium is a routine chore for the particle cosmologist. The diff erential equa- 
tion governing the abundance of a species follows from the Boltzmann equation, 
and depends upon two pieces of input physics: the expansion rate as a function 
of temperature and the annihilation rate as a function of temperature. 2 Once the 
particle species and its interactions are specified the annihilation rate is precisely 
determined in terms of the number density of the species and the temperature of 
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the Universe. The expansion rate as a function of temperature is another mat- 
ter. In the standard, radiation-dominated FRW cosmology, the expansion rate 
at early times (t < tEQ ) is given by 


H 2 


2 

R\ _ 87 xG _ 47T 3 g»r 4 
R ) 3 P ~ 45 m|, ’ 


( 1 ) 


where R is the scale factor of the Universe, p is the total energy (for a thermal 
bath of relativistic particles p T = £r,7r 2 T 4 /30), and g * counts the effective number 
of ultra- relativistic degrees of freedom (1 for each relativistic bosonic degree of 
freedom and 7/8 for each relativistic fermionic degree of freedom). Having made 
this assumption, the path to determining the relic abundance — usually expressed 
as the ratio of the number density n of the species to the entropy density s = 
2g^ 2 T 3 /45 — is a tried and true one. 

The crucial uncertainty in determining the relic abundance is the assump- 
tion that the Universe is radiation dominated at freeze out. The concordance 
of the predictions of primordial nucleosynthesis with the observed light-element 
abundances provides strong evidence that the Universe was indeed radiation dom- 
inated at an age of about 1 sec. However, there is little or no evidence that 
requires the Universe at times earlier than about 1 sec to be radiation dominated. 
Moreover, most of the interesting thermal relics would have decoupled at such 
early times. 

Given the importance of relic particles both to cosmology and to particle 
physics — particle relics may provide the bulk of the mass density of the Universe 
thereby explaining the nature of the dark matter, and they might provide the 
first evidence for new physics beyond the standard model of particle physics 
we decided to study in detail three nonstandard — but plausible — possibilities for 
the expansion rate of the Universe around the time of freeze out. Our first 
example involves the geometry of spacetime: If the expansion of the Universe is 
not isotropic, then the volume-expansion rate at fixed temperature exceeds that 
for the standard case and freeze out occurs at a higher temperature leading to 
a larger relic abundance. Here we explore a particularly simple and interesting 
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example: a Bianchi I model where the effects of the anisotropy on the volume- 
expansion rate can be quantified in terms of an anisotropy-energy density that 
decreases as R ~ 6 . For this model (and other similar models) the anisotropy 
simply decays without leaving a trace, and the only lasting effect is to enhance 
the abundance of the thermal relic. 

In the second example we consider, at early times the energy density is dom- 
inated by a massive, nonrelativistic particle species. Again, the expansion rate 
for fixed temperature is increased, leading to an earlier freeze out at a larger relic 
abundance. Of course there is every evidence that the Universe only became 
matter dominated relatively recently and so the nonrelativistic particles would 
eventually have to decay, producing entropy and diluting the relic abundance 
of the relic. As we shall show the net effect is to decrease the relic abundance. 
(We note that this possibility is different from the one involving particle decays 
in which the only effect of the decaying species is to produce entropy, in which 
case the relic abundance is decreased precisely by the amount of the increase in 
entropy. 5 ) 

In the third example we use the Brans-Dicke-Jordan theory of gravity instead 
of general relativity. Here the analog of the Friedmann equation contains a kinetic 
energy term for the Brans-Dicke scalar scalar field which decreases as i? -6 and 
of course increases the expansion rate for fixed temperature. As in the case of 
anisotropic expansion, the only lasting effect is to enhance the abundance of the 
thermal relic. 

The motivation of this work then is to assess the reliability of the standard 
estimate for the relic abundance of a stable particle species that was once in 
thermal equilibrium by considering three nonstandard possibilities for the evo- 
lution of the Universe at early times (t < 1 sec). The outline of our paper is 
as follows: In Section II we briefly review the formalism for calculating the relic 
abundance of a species and the standard result. In the next three Sections we 
consider the nonstandard possibilities mentioned above and how they affect the 
relic abundance of a stable particle species. In the final Section we put our work 
in perspective with some concluding remarks. 
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II. REVIEW OF THE STANDARD RESULT 


To obtain quantitative results for the relic abundance of a stable particle 

species X (and its antiparticle X ) one solves the Boltzmann equation that governs 

2 

the number density of the species: 


— + 3 Hn = -{av)(n 2 - n| Q ). (2) 

Here n(t ) is the number density of species X at time t, tieq is the equilibrium 
number density at time t when the temperature of the plasma is T, ( crv ) is the 
thermally and spin-averaged cross section times relative velocity for A A anni- 
hilation, and H is the Hubble parameter. We assume that there is no particle- 
antiparticle asymmetry so that the number density of antiparticles is also n. (It 
is easy to extend this formalism to apply to the case where there is a particle- 
antiparticle asymmetry; see Griest and Seckel. ) Exact solutions of this particular 
example of the Ricatti equation do not exist; however, an approximate [accurate 
to 0(5%)] analytical solution is easily obtained. Since we will follow this ap- 
proach in the nonstandard models we will briefly review that solution here. 

In the absence of entropy production the entropy per comoving volume (5 = 
R 3 s) is constant, and we use the entropy density 


Z7r rp 3 

5 = ls S - T ’ 


as a fiducial quantity and rewrite the Boltzmann equation in terms of Y = n/s, 
which corresponds to the number of X particles per comoving volume. [Actually, 
the g * in Eq. (1) is in principle different from that in Eq. (3); however, in practice 
they are very nearly equal at temperatures greater than an MeV. See Ref. (3), 
pp. 65-70.] Furthermore, since the quantities of interest depend explicitly on 
temperature rather than time, we use the quantity x = m/T, instead of time as 
the dependent variable. Doing so, Eq. (2) becomes 


dY 

dx 


< (crv)s 


( y 2 - y | Q ), 


( 4 ) 


where of course (cv), H , s , and Yeq are all functions of x . (Actually, there is 

6 8 

another term involving the derivative of #*; 5 however, this term is small and 
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to a good approximation we can fix g+ at its value around freeze out.) Here 
Ieq = ngg/s and in the nonrelativistic limit (x 3) is given by 


45 /ic\ 1/2 g 


Y EQ(x) = — 4 {-) 


—x 3 /~e~ x = 0.145— x 3/,2 e -1 . 


(5) 

9* 9* 

In the case of interest, at freeze out XX particles are moving at nonrelativistic 
velocities and the cross section is proportional to v 2n (n = 0 corresponds to s- 
wave annihilation, n = 1 corresponds to p-wave annihilation, etc.) so we can 
write ( av } = crox -n . Furthermore, s oc x -3 and H oc x -2 , cf. Eq. (1), so the 
Boltzmann equation becomes 




A 


r 2 


where we have defined 

A = 


{ [crv)s 

H 


x 


x = l 


= Q.264g¥~mp\mxcro. 


( 6 ) 


(7) 


To solve Eq. (2) we follow Ref. 7 and consider the differential equation for 
A = Y — Teq> the departure from equilibrium, 

' A ' = ->eq - (2Y eq + A) , (8) 

where the prime denotes d/dx. At early times when the annihilation rate Tann 
is much greater than the expansion rate H (x <C x/), the X abundance tracks 
its equilibrium value very closely so that A <C Yeq and A' <C Yg q, and 

_n+ 2yt 

A ~ - ( 9 ) 

~ A(2Ye Q + A) ' ( • ) 

At late times (x >■ x _/■), Y tracks Yeq very poorly; therefore A~7> Yeq and 
Y^q <A', so that 

a ' = -^F2 a2 - (10) 

7 

Upon integrating Eq. (10) from x = x j to x = oo, we obtain 

Too = Aoo = — P*/ +1 + 0(x n f), (11) 

where Xf is determined from Eq. (9) by A(xy) = cYeqC^/) and c is a numerical 
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constant of order unity that serves to define the epoch of transition between the 
regimes mentioned above. Solving for xj iteratively, the result is 

Xf ~ In [(2 + c)Xac] - (n + ^ In {In [(2 + c)Aac]} , (12) 

where a = 0.145 (g/g*). Note that the final abundance only depends logarthically 
upon the value of c. The best agreement between the analytic result and a 
numerical integration of Eq. (4) obtains for c(c + 2) = n + 1. To obtain the 
present mass density contributed by the relic, Yoo is multiplied by the mass of 
the relic and the present entropy density so = 2970 T| 75 cm -3 , where T2.75 is the 
CMBR temperature in units of 2.75 K. 


III. FREEZE OUT IN A SHEAR-DOMINATED UNIVERSE 


As stated earlier, there is little or no evidence that the Universe before the 
time of big bang nucleosynthesis had to be homogeneous or isotropic, although 
this is the standard assumption and is certainly well motivated. The simplest 
nonstandard cosmological models are homogeneous but anisotropic; these are the 
Bianchi (and Kantowski- Sachs) models which are classified according to the Lie 
algebras of their isometries. 9 The metric of the simplest of these models, the 
Bianchi type I spacetime, is 

ds 2 = _ dt 2 + Rltfidx 1 ) 2 + R\{t){dx 2 ) 2 + Rl(t)(dx 3 ) 2 , (13) 


where the Ri are the scale factors of the three principal axes of the Universe. The 
Einstein equations for this metric lead to the analog of the Friedmann equation 
for the volume expansion rate H of the Universe: 


H 2 


1 

9 
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3m pj 


(pr + Ps), 


(14) 


where we assume that the matter content is the usual thermal bath of radiation 
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at temperature T and the shear “energy density” is defined to be 

P, = jjjj [(#1 - H 2 ? + (Hi - H,? + (Ht- Hi?} . (15) 

Here V = R 1 R 2 R 3 is the “volume scale factor,” R = U 1 / 3 is the mean-scale factor, 
and the Hi = (R t /R,) (no sum) are the expansion rates of the three principal 
directions. As is manifest from Eq. (15), the shear-energy density is proportional 
to the amount of anisotropic expansion. Note that since we axe always free to 
relabel our comoving coordinates, differences between the various R x 's have no 
physical meaning; only differences in the expansion rates axe meaningful. From 
Eq. (15) it also follows that |i?i| < 3 H (we use absolute-value signs because 
at very eaxly times, when p T is negligible, the spacetime becomes the Kasner 
spacetime in which one of the spatial dimensions must be contracting). 

In general, the red shift suffered by a particle as the Universe expands will 
be direction dependent. For example, for a particle moving in the x-direction, 
p oc RJ 1 . Provided that the interaction rate of the thermal bath of particles is 
much larger than H, particle distributions will remain isotropic and the mean 
m ome nta will red shift as R~ l . In this case, the remaining Einstein equations 
become (for i ^ j): 

j t \n\H i -H J \ = -3H = -3j t (lnR), (16) 

which implies that p s oc Rr\ Therefore, the shear-energy density falls off faster 
than the radiation-energy density and the anisotropy in a Bianchi I Universe 
simply decays without leaving a trace. 10 

For the freeze-out calculation we are interested only in the expansion rate 
H and not the detailed form of the anisotropy as given by the Hp, therefore, we 
use the fact that p s oc R ~ 6 and constancy of the entropy per comoving volume 
(g t R 3 T 3 = const) to express the shear-energy density in terms of the plasma 
temperature T. We define the temperature at which p T = p s to be T e . For 
T > T e the Universe is shear-dominated, H oc R~ 3 and R oc f 1 / 3 ; for T<T e 
the Universe is radiation-dominated, H oc R ~ 2 and R oc t 1 ! 2 . The temperature 
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T e quantifies the size of the anisotropy energy: Smaller values of T e correspond 
to larger anisotropy energy density at fixed temperature. We then write the 
shear-energy density in terms of the radiation-energy density, 

= *(T.) (^) =^(^) • (17) 


where g% is the value of g * at T e . 

In order to avoid interfering with the successful predictions of big bang nucle- 
osynthesis, we must be sure that the shear-energy density is sufficiently small at 
the time of primordial nucleosynthesis. The shear contribution to the energy den- 
sity would speed up the expansion rate thereby increasing the 4 He production. 
Assuming that the only contribution to the energy density comes from radiation, 
concordance of the outcome of nucleosynthesis with the observed abundance of 
4 He requires that g*(T ~ MeV) < 12.5.^ To assess the effect of the shear-energy 
density, we write the total energy density as 


P = 


—q ef f 
3(r* ’ 


(IS) 


where 

(i9) 

the primordial nucleosynthesis constraint is then g e J^ < 12.5. (In terms of the 
number of light neutrino species, this is equivalent to N v < 4. Recent SLC 
and LEP results have confirmed this constraint, determining that the number 
of neutrino species lighter than about 40 GeV is 3.2 db 0.2. ) If the r neutrino 

is light (m„ r < few MeV), we know that g+ is at least 10.75. Requiring that 
< 12.5 then leads to the constraint: T e > 2.5 MeV. 

At this point, we can see how shear can increase the relic abundance of a parti- 
cle species. In Fig. 1 we plot the expansion rate H in a shear-dominated Universe 
as a function of temperature T. At low temperatures (late times), H ~ T 2 , while 
at high temperatures (early times), H ~ T 3 . The broken curve shows the ex- 
pansion rate H si d with no shear. The equilibrium number density of X is 
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proportional to T 3 at high temperatures (T >• m) and falls exponentially at low 
temperatures (T <C m). For the case that the thermally and spin- averaged cross 
section times relative velocity, (crv), for XX annihilation is constant, the annihi- 
lation rate, I"ann = n EQ,{ av ) n EQi is also shown in Fig. 1. Roughly speaking, 
at the freeze-out temperature Tf, defined by Tann = H, annihilations freeze out, 
and the number of X’s per comoving volume “freezes in,” at approximately its 
value at Tf. If Tf > T e , the expansion rate in the shear-dominated Universe 
is much greater than that in the standard radiation-dominated model, and the 
annihilations freeze out earlier when the abundance is greater. Since Tann de- 
creases exponentially around freeze out, the freeze-out temperature for the two 
cases (shear and standard) is nearly the same. Moreover, because tieq oc Fann? 
the relic abundance in a shear-dominated model {Tf > T e ) is enhanced roughly 
by a factor H(T f )/H sti (T f ) ~ T f /T e . 


To obtain more quantitative results for the relic abundance, we must solve 
the Boltzmann equation using the expression for H which includes the effects of 
anisotropic shear. In the Appendix we show that the Boltzmann equation used 
in the standard cosmological model, Eq. (2), is also valid in the Bianchi I model 
as well. For an anisotropic-universe model the expansion rate is conveniently 
written as H = H s td(z 2 + x 2 ) 1 ' 2 j x where U s td is the standard model expansion 
rate and 


m 




( 20 ) 


Doing so, the Boltzmann equation, cf. Eq. (4), becomes 


dY_ 

dx 


A 

X n+l {x1 + x 2 ) 1 / 2 


(r 2 - Ieq)- 


( 21 ) 


The differential equation for A then becomes 


A' = - 


As before, at early times, 


A 


x 


n+1 


(** + X?)>/2 A ^ + A ) ' 


( 22 ) 


x“ + V +I 2) i/sy^ 
■1(2V"eq + A) 


(23) 
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and at late times 


A 


(24) 


A' = 


X n+ 1 (x 2 + X 2 ) 1 / 2 


For a model with no shear (x e = 0) we recover the results of the previous 
Section. For a model universe where shear is important at freeze out (x e > x/), 
an exact closed form solution for Eq. (24) for arbitrary n is not simple to write 
down. However, if x e >• Xf, a good approximation for Y+, may be obtained by 
integrating in the interval Xf < x < x e assuming H = H s td/x and in the interval 
x e < x < oo assuming H = i? s td- Doing so we obtain 


Yoo 


nx e x n j 

A 


1 + 0 



(25) 


for n^O, and 


Xe 

Aln(2 x e /xfY 


(26) 


for n = 0. Assuming that freeze out occurs while the Universe is shear dominated 
(xf <C x e ), the equation for Xf is given by 


Xf ~ In [(2 + c)acAx e l ] 


We see that x f decreases roughly by only an additive factor of ln(x e / x f) justifying 
our previous assertion that the freeze-out temperature is nearly the same in a 
shear-dominated or radiation-dominated model. 


Defining an enhancement factor 


£ 


(tV 12 Tl 
WJ T.’ 


(28) 


we see that if the particle-antiparticle annihilation is primarily s-wave, the relic 
abundance in a shear-dominated Universe is increased roughly by a factor £/ ln£ 
over that in the standard case. If the annihilation is primarily p-wave, which is 
often the case for Majorana particles, the enhancement is roughly 0.5 £. This 
result is particularly interesting for Majorana particles (e.^., Majorana neutrinos, 
photinos, higgsinos, etc.) since an enhancement in the relic abundance due to a 
particle-antiparticle asymmetry is not possible for self-conjugate particles. 
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As an example of current interest, we may apply our results to a Majorana 

neutrino of mass 60 GeV. Since such a neutrino is heavier than half the mass of 

the Z°, the decay — * vv is kinematically forbidden. Thus, such a fourth gen- 

13 

eration neutrino is not exluded by the recent SLC/LEP results. Furthermore, 
since it has only “spin-dependent” couplings to nuclei, its elastic scattering cross 
section is too small for it to be ruled out by the results of germanium ionization 
experiments. 15 Even so, it is not generally considered a candidate for the pri- 
mary component of the dark matter in the galactic halo since its abundance — as 

determined by standard calculations — is small (f2„ 1) and cannot be increased 

1 6 

by introducing a particle-antiparticle asymmetry. 

In Fig. 2, we show the results of a numerical integration for Y as a function 
of x for the cases of x e = 0 (standard model), x e = 1300 (shear-dominated model 
with T e ~ 120 MeV) and x e = 13000 (T e ~ 12 MeV). For all three cases, the 
numerical results agree with the analytic results, Eqs. (25) and (27), to within 
5% [using c(c + 2) = n + 1]. In the standard cosmology, the resulting value of 
Q, u h 2 is 1.1 x 10 -3 — too small for heavy neutrinos to be the primary component 
of the galactic halo. We find that for T e ~ 120 MeV the present mass density is 
increased to Q,„h 2 ~ 0.021, a value comparable to that known to be contributed 
by the halos of spiral galaxies; and for T e ~ 12 MeV the present mass density is 
f l„h 2 ~ 0.19 which is about right to close the Universe. 

From this we conclude that some stable particle species that have not been 
considered dark matter candidates due to their small relic abundances could 
indeed still be dark matter candidates. Only dark matter search experiments, 
such as germanium ionization experiments or future bolometric detectors, can 
definitely rule out a particle species as being the primary component of the 
galactic halo. Moreover, if an “unlikely” particle relic is discovered, cosmologists 
would have to significantly alter their current notions of the first second of the 
Universe’s history. 

Finally, we mention a possibility suggested by Misner and others, the “decay” 
of anisotropy into radiation due to dissipative processes. 10 If this occurs, then 
the anisotropy-energy density at freeze out could have been much larger than the 
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upper bound imposed from nucleosynthesis (T e > 2.5 MeV), provided that the 
dissipation took place before the epoch of nucleosynthesis. Naively one might 
expect that the enhancement of the relic abundance could be arbitrarily large; 
however, this is not correct. The entropy produced by the dissipation of the 
anisotropy will dilute the relic abundance. Ignoring factors of order unity, if the 
anisotropy is dissipated at a temperature Tp, the ratio of entropy per comoving 
volume after dissipation to that before dissipation is (Tp/T e ) 3 / 2 ; the relic abun- 
dance is reduced by this factor. The net enhancement over the standard result 
is about 


(IlY 2 it 

\T d J T d ' 


(29) 


Since T e < Tp < Tf, this factor can be greater than one; however, since 
T d > O(MeV), the enhancement can never be as great as the maximum en- 
hancement allowed by our previous analysis where there was no dissipation and 
T e > O(MeV). 


IV. FREEZE OUT IN A MATTER-DOMINATED UNIVERSE 

Next, consider a model where the energy density of the Universe at the freeze 
out of particle species X is dominated by some other massive particle species 0 
which subsequently decays. Before the decay of 0, the Friedmann equation is 

jr, -5sfc( T ‘ + " 34 }' (30) 

where M is a very large (M » Tj) parameter with dimensions of mass: Specifi- 
cally, M = 4m©n©/3s where m© is the mass and n© is the number density of 0 
particles. Once again, since freeze out occurs roughly when TanN = H and the 
number density is proportional to TanNj the relic abundance is apparently en- 
hanced by roughly (M/Tf) 1 / 2 . To be more precise, since M > Tf the expansion 
rate in this model is if = ff st d(Mr/m) 1/2 where if st d is the expansion rate in 
the standard cosmology; therefore, the Boltzmann equation for this case is given 
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by Eq. (4) with the substitutions A — ► A (m/M) 1 / 2 and n — > n + 1/2. Making 
these substitutions we can “read off” the solution from Eq. (11); the result is 


Yoo = 


(n + 3/2 )M n+3/2 

— £ _ 

Am f 


+ 0(x n + 1/2 ) 


(31) 


which agrees with our rough guess of the enhancement factor. 

However, this is not the whole story. The subsequent decays of 0 particles 
occur out of equilibrium and produce a large amount of entropy thereby lessening 
the previous enhancement. In fact the net result is a reduction in the relic abun- 
dance relative to the standard case. To see this, suppose that the temperature 
at which the © particles decay is To (which of course is less than T/), then the 
ratio of entropy per comoving volume after decay to that before decay is roughly 
(M/Tz?) 3 / 4 ( see Ref. (3)); therefore, the final relic abundance is roughly a factor 

T 3 V /4 

< 32 > 

times that in the standard case, given by Eq. (11). Since Tp < Tf < M, no 
enhancement in the relic abundance is possible; rather, the relic abundance is 
reduced. We could attempt to circumvent the entropy-production problem by 
supposing that ©-particles decay into some non-interacting, inert species that 
does not contribute to the “visible” entropy density. However, in any interesting 
case, the additional relativistic degrees of freedom would exceed those allowed by 
primordial nucleosynthesis. 


V. BRANS-DICKE-JORDAN COSMOLOGY 

There has been renewed interest interest in alternative theories of gravity, 
particularly those in which the gravitational “constant” varies as it is does in the 
Brans-Dicke-Jordan theory. Much of this interest owes to the advent of extended 
inflation, a variant of old inflation in which the “graceful exit” problem is 
solved. Although it now appears that extended inflation in the Brans-Dicke- 
Jordan theory is not viable, as the isotropy of the microwave background 
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requires the Brans-Dicke parameter uj to be less than about 30, while solar- 

20 

system experiments require that u ; > 500, variants of the Brans-Dicke- Jordan 
theory may still lead to successful inflationary scenarios. In this section we will 
show that a cosmological model based on the B rans- D icke- Jordan theory with 
lj > 500 allows for significant enhancement in the abundance of a thermal relic. 
Since many of the scalar-tensor theories currently under consideration resemble 
Brans-Dicke-Jordan theory (with a variable uj), we expect that our results may 
generalize to these theories as well. 

The Brans-Dicke-Jordan theory of gravitation is the scalar-tensor theory 
that can be derived from the action 



— + uj 


$ 


~b 167r.C ma tter > 


(33) 


where 7 Z is the curvature scalar, and the real scalar field $ has dimensions of 

mass squared and sets the value of the gravitational constant: G = 4> _1 ; for 

this reason, $ must necessarily be greater than zero. Since $ is a dynamical 

field, one expects the gravitational constant to evolve with time. The quantity 

uj is the dimensionless Brans-Dicke parameter; in the limit that uj — ► oo, the 

scalar-tensor theory reduces to general relativity. While the scalar-tensor theory 

becomes much less attractive for uj » 1 , it still provides a simple example of the 

kind of different gravitation theory that might arise as the low-energy limit of 
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some superstring models. 

Specializing to the Robertson- Walker line element and for simplicity to a 
spatially-flat model, the equations of motion for the scale factor R(t) and for $ 
axe: 

!(„*’) = - 4 * 3 ; (34 ) 



87T 

3 + 2uj 


(p- 2p)R 3 ] 


(35) 


H 2 



8t rp uj 

"3$" + ? 




(36) 
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$ 

ff = '2? + 


(2u + 3) / $ \ Stt/> 
3 2$ / + 3$" 


1/2 


( 37 ) 


where as usual /9 is the energy density (of all the fields other than $) and p the 
isotropic pressure. In going from Eq. (36) to Eq. (37), we have considered only 
the positive root, as we are interested in expanding Universe models. 

Note that in Brans-Dicke-Jordan cosmology there are two additional bound- 
ary conditions that must be specified: the values of $ and $ at some epoch. 
Since the theory must closely resemble general relativity today, the present value 
of $ must be equal to G — 1 (for large ui; see Weinberg, Ref. 18): $o — . 

That effectively specifies one of the boundary conditions. The other, involving 
the value of at some epoch, still remains to be specified. 

The Brans-Dicke-Jordan analogue of the Friedmann equation, cf. Eq. (37), 
differs from the usual one in two regards: First, the gravitational constant is 
given by $ -1 ; second, there is an additional contribution to the energy density 
that involves the kinetic energy of the $ field. It will be useful to consider the 
the ratio of the ^-kinetic term to the usual energy density term: 


2 

_2w + 3($\ /Si T P (2u> + 3)$ 2 

r=_ 3“(2fj /sF = 32rf>p ; (38) 

as we shall see the ratio r decreases with time: r oc R~ 2 (when p is radiation 
dominated); and r a const /(Inf) 2 (when p is matter dominated). Having defined 
r we can rewrite Eq. (37) in a very suggestive form: 


H = 



(39) 


where the upper sign applies for > 0 and lower sign for <£ < 0. In Eq. (39) 
the two modifications to the usual Friedmann equation are manifest: For r / 0 
the presence of the $ field speeds up the expansion rate; and if ^ G~ l the 
expansion rate is also changed. 
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A. Energy density dominated by relativistic particles 

To begin, let us consider the case where the energy density of the Universe 
is dominated by relativistic particles, which is what one expects at very early 
times. In this case, p — pj 3, so that 

^-($fl 3 ) = 0 =* $R 3 = B\ 

at 

I W*)— 55* * pRi = A '’ 

where A and B are numerical constants. (We will neglect the slight variation of 
A that occurs because g* evolves.) In terms of A and B , r is given by 

(2w + 3 )B 2 
r ~ 32tt AST? 2 ‘ 

We see that the $ boundary condition can be set by specifying the value of B, 
or equivalently the value of r, at some epoch. During the radiation-dominated 
epoch the value of $ does not change very much, so that r oc R~ 2 ‘, stated another 
way, the additional energy density associated with the $ field red shifts as R~ & . 

It is simple to integrate the equations of motion, Eqs. (35) and (37), to obtain 
$ as a function of R : 24 



$(i?) = c 


(1 + r 1 ) 1 / 2 — 1 
(1 + r 1 ) 1 / 2 + 1 


-| ±-y/3/(2a/+3) 


(41) 


where the upper sign applies if B > 0 and the lower sign if B <0, and the 
constant of integration C manifests the freedom one has to specify the value of 
$ at some epoch. At early times, corresponding to small R and large r, 

(4r) T v /3 /( 2 “+ 3 ) . 


while at late times, corresponding to large R and small r, 

$ -> C( 1 - 2Vr) ± v' 3 /( 2w + 3 ). 


Since u 1, at early times when r 1, the value of $ slowly increases for B > 0 > 
(decreases for B < 0); once r ~ 1, the value of $ asymptotes to the value $ = C 
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(regardless of the sign of B). When r <C 1 — dynamics controlled by the energy 
density in radiation — $ ~ C, and r a R~ 2 . And of course, the scale factor of 
the Universe grows as f 1 / 2 . In this regime the expansion of the Universe behaves 
as is if there is an additional form of energy density that decreases as R ~ 6 — just 
like shear. 

In order that the successful predictions of primordial nucleosynthesis not be 

12 

upset, r must be less than about 0.2 when t ~ 1 sec and T ~ 1 MeV; this 
constrains the initial value of 4>. (Moreover, we must also insure that the value 
of $ does not differ from its present value by more than about 20%; as we shall 
see below, this only requires that > 50.) The constraint is 

£> 0.2 ( 2^ 3 ) h BBH ~ 2 ^^+ 3) ° r ~ ^bbnI*bbh)I \/2(2kJ + 3). 

(42) 


B. ^-dominated expansion dynamics 

Since r evolves as R ~ 2 , at early times the dynamics of the expansion will 
necessarily be dominated by the $ field. For r >> 1, the equation for the expansion 
rate of the Universe becomes, 


H 


4 


2u> 3 


1/2 


Tl 


$ $ 


(43) 


where (3 = [\J(2u -f 3)/3 l]/2 > 0, and the upper sign applies if 4 > 0, 

while the lower sign applies if $ < 0. Assuming that the energy density of the 
Universe is still dominated by relativistic particles, this equation is supplemented 
by $ = B/R 3 . 

These equations are straightforward to solve: 

R(X$ ±f} i? oc $ oc *±i/(30±i) ; 



R— 6=T=2//J 


T CK 




Since j3 ~ y / o>/6 is expected to be large (greater than about 10), it follows that 
R oc t 1 / 3 , $ 2 /<$ 2 oc i? -6 , r oc i? -2 , and $ increases slowly with time for $ > 0 
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(decreases for $ < 0). That is, during the ^-dominated phase the Universe 
behaves like a FRW model whose expansion dynamics axe controlled by a form 
of energy density that decreases as R ~° — just as in a shear-dominated model. 


C. Energy density dominated by nonrelativistic matter 

At an age of about t£Q ~ 4 x 10 10 sec and temperature of about Teq ~ 10 
eV, the Universe becomes matter dominated. Based upon the nucleosynthesis 
bound, we can infer that teq ~ 10 ~ iq tbbn ~ 10 -11 and 


$EQ = C 


1 T 2 


2 u) -p 3 


1/2 


T'EQ 


1 $■£<? ~ £1 < 10 11 
C ~ y/Z ' 


To a very good approximation, the Universe will behave as an ordinary matter- 
dominated FRW model and R/Req = (t/tEQ) 2 ^- During the matter-dominated 
epoch pR? = const, and it is convenient to express the value of that constant as 
PeqR'eq- Thus, the evolution of $ is given by 

= iSs'**** (44) 


L ,2 4 $EQtEQ] 1 4 $EQ 1. 

m - [*«»*£* - 3(2u , + 3) j + 3 2ZT3T 


*( () = + 3 & + 3 ) ln(</<£0) + {**>**, ” 3 ( 2 ^ 3 ) 


(45) 

1 - ■ (46) 


From Eq. (46) we can find the value of $ at the present epoch (t = to ~ 10 7 f£Q): 


$0 = $EQ + 3 ^2u + 3) [ ln (W*£Q) “ !] + $EQtEQ-, (47) 

from our constraint to $ eq it is simple to show that the term involving $ eq is 
negligible: ^EQ^EQ ^ 10 -4 $£q/2(2u> + 3) 1 / 2 . In order that $0 not differ from 
$EQ by more than about 20%, u> must be greater than about 50, which is not as 
stringent a bound as that provided by the solar system experiments. Finally, it is 
simple to see from Eq. (45) that in the matter-dominated epoch r oc const/(lnf) 2 , 
while $ grows logarithmically with time. 
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To summarize, primordial nucleosynthesis constrains r to be less than about 

0.2 at the epoch of nucleosynthesis and u to be greater than about 50. The 

constraint to r provides information about the the initial value of $. At very 

early times the dynamics of the Universe are necessarily controlled by $ since 

r oc R~ 2 . The transition to ^-dominated expansion dynamics will occur at 

a temperature of about ~ r BBN MeV, which could be as low as 3 MeV. 

During the phase when $ controls the dynamics of the expansion, the Universe 

behaves like an ordinary FRW model whose energy density is dominated by a 

form of energy that decreases as i? -6 : R oc i 1 / 3 . This has implications for the 

relic abundance of a thermal relic that freezes out at a temperature greater than 

about 3 MeV, which will address below, as well as for coherent axion production 
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and for baryogenesis, which we will address elsewhere. 

The analysis of the “freeze in” of the relic abundance of a stable particle 
species that freezes out at a temperature Tf > is identical to that in the 
previously discussed shear-dominated model. That is, the relic abundance is 
increased, relative to the standard case, by a factor of £/ln£ for s-wave, or 0.5<f 
for p-wave, where 

£ = (9J9? ) 1/2 ^; 

where gf is the value of g* when r = 1. We should point out that for the 
simplest Brans-Dicke-Jordan extended-inflationary model, the Universe enters 
the radiation-dominated epoch directly at the end of inflation bypassing a d>- 
dominated epoch so no enhancement in the abundance of a thermal relic can 
occur. 24 However, the details of extended inflation axe fax from being com- 
pletely understood— including whether or not inflation took place — and so a 
^-dominated epoch is an interesting cosmological possibility. 

As is clear from this Section and the previous two, what is required to enhance 
the relic abundance of a particle species is that the Universe at early times be 
dominated by a form of energy density that decreases faster than R ~ 4 . In this 
case, this component of the energy density can dominate the energy density 
at freeze out and then conveniently disappear before primordial nucleosynthesis 
without leaving a trace. Shear in a Bianchi I model and the kinetic energy of 
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the Brans-Dicke field provide two examples where the additional energy density 
is proportional to R~ . 

There are more. The energy density of any fluid for which the pressure 
exceeds the energy density divided by 3 will decrease faster than i? - . The 
extreme case is p = p, which used to be discussed as an equation of state for 
the Universe at very early times. A homogeneous scalar field with Lagrangian 
density 

behaves like a perfect fluid with equation of state p = (n — 2)p/(n + 2) as it 
oscillates about the minimum of its potential (|</>| = 0). In so doing the associated 
energy density p oc i?~ 6n /( n+2 ). 27 For n > 4, the energy density of such a field 
decreases faster than il -4 . For {3 = 0 or n — * oo, <f> is a massless, free scalar field 
and p cx R~ 6 . (An interesting example of a massless scalar field is a Goldstone 
mode. Consider a complex scalar field (f> = ipexp(i9) with a “Mexican-hat” 
potential, where ip and 6 are real scalar fields. Suppose that the magnitude of ip 
is fixed by spontaneous symmetry breaking; the phase 6 is a massless Goldstone 
mode, which can spin around the brim of the hat. In fact, this is precisely what 
occurs in a recently suggested scenario for baryogenesis. ) 

Along similar lines as the Brans-Dicke- Jordan theory is the possibility that 
the gravitational constant varies because it is related to the size of some extra, 
compactified dimensions. In many Kaluza-Klein and superstring theories the 
gravitational constant G varies as G = G to< j ay (T/i t oday) jD where L is the scale 
factor of D compactified dimensions. If one assumes that all 3 + D spatial 
dimensions were of comparable magnitude at early times, it suggests that early on 
the gravitational constant was larger than it is today. The strongest constraint to 
the variation of G is that imposed by primordial nucleosynthesis, which implies 
that by the epoch of nucleosynthesis the value of G differed from that today by 
less than about 20%. However, there are no stringent constraints to the value of 
G at earlier times. If it were very different than its present value — and larger — 
then the expansion could have been faster than in the standard cosmology. Since 
Foo oc rripj oc G 1 / 2 , cf. Eqs. (7) and (11), we would expect the relic abundance 
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VI. CONCLUDING REMARKS 


to be increased. 


While much of the activity in cosmology these days involves the study of 
the earliest moments of the Universe, we have precious few probes of those early 
times. A class of potential probes are thermal particle relics — stable particle 
species that were once in thermal equilibrium. Already such relics have received 
a great deal of attention, particularly as candidates for the dark matter. Ther- 
mal relic dark matter candidates include heavy neutrinos, neutralinos, and light 
neutrinos to mention three of the most interesting possibilities. It goes with- 
out saying that the discovery of such a relic would be of enormous importance 
to cosmology; in addition, the discovery of any of the aforementioned particle 
species would be of equal importance to particle physics, providing evidence for 
new physics beyond the standard model of particle physics. 

The calculation of the relic abundance of a particle species has become a very 
routine task for the particle cosmologist. In this paper we have addressed the 
crucial and untested assumption in the calculation: the temperature dependence 
of the expansion rate of the Universe. In nonstandard cosmological models where 
the energy density of the early Universe is dominated by nonrelativistic matter, 
anisotropy, or the kinetic energy of a scalar field, the relic abundance can be 
significantly different. In the case of the energy density being dominated by non- 
relativistic matter the relic abundance is ultimately smaller than in the standard 
case, due to the entropy produced by the eventual decays of the nonrelativistic 
particles. In the case of a Universe that is shear-dominated, or ^-dominated, 
early on, the relic abundance can be greatly enhanced owing to the fact that 
the expansion rate for a given temperature is larger, which leads to a freeze out 
at a higher temperature and a larger abundance. We remind the reader that in 
spite of the fact that the standard, radiation- dominated FRW model is very well 
motivated, there is no direct evidence that excludes the possibilities that we have 
discussed here. (We do mention that the levels of shear that are interesting for 
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3 1 

our purposes are definitely incompatiable with the inflationary Universe. ) 


The fact that the relic abundance of a particle species can be greater than 
the canonical abundance is of no small interest to those involved in dark mat- 
ter searches. It is well known that the results of the standard relic abun- 
dance calculation can be decreased by phenomena — such as inflation, an elec- 
troweak or quark/hadron phase transition, or out-of-equilibrium decay of a mas- 
sive particle— that produce a significant amount of entropy after freeze out. If 
the canonical calculation indicates that the relic abundance is too small for the 
species to be the primary component of the galactic halo, any of these entropy- 
producing processes only make the conclusion that much stronger. However, if 
freeze out occurs in a shear- dominated or ^-dominated epoch, the relic abun- 
dance is enhanced. Thus, some particle species that are not interesting dark 
matter candidates according to the standard calculations may indeed be inter- 
esting dark matter candidates. Perhaps one should take the empirical view that 
a particle dark-matter candidate should only be ruled out by null results in dark- 
matter searches or accelerator searches. Our work also implies that the rates for 
indirect signatures, such as high-energy neutrinos from particle dark matter anni- 
hilations in the Sun or Earth, or positron-line or 7-ray line radiation from particle 
dark matter annihilation in the halo, could be significanty larger than expected. 
Finally, the discovery of one of these “unlikely” particle relics in the galactic halo 
would force us to reconsider our current view of a radiation-dominated Universe 
at times earlier than about 1 sec. 
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APPENDIX 


Here we show that the Boltzmann equation, cf. Eq. (2), used to describe the 
evolution of the number density of a species in the FRW model is also valid for 
the Bianchi I model considered in this paper. To do so we follow the discussion 
in Ref. (3). 

The evolution of a particle’s phase-space distribution /(p /1 , .t / 1 ) is governed 
by the Boltzmann equation, which can be written as 


£[/] = c[/], 


(4S) 


where C is the collision operator and L is the Liouville operator 


CULJLVk lO 




(49) 


where are the usual Christoffel symbols. For the Bianchi models, the phase- 
space density is spatially homogeneous so / is a function of time t = x° but not 
space. A crucial assumption is that scattering interactions are occuring rapidly 
enough so that the particle species remains in kinetic equilibrium. Provided 
that this is the case, the phase-space density / is isotropic and only depends 
on the magnitude of the momentum, or equivalently the energy E. With this 
assumption, for the Bianchi I metric, cf. Eq. (13), the Boltzmann equation 
becomes 


L[/(E,<)] = E %- [^(P 1 ) 2 + «2J?2(p 2 ) 2 + RsJMp 3 ) 2 ] J|. (50) 


Since the number density of the species is 


n(t) = J d 3 pf(IS,t), (51) 

the equation for the evolution of the number density is obtained by multiplying 
Eq. (48) by gd?p/(2Tr)E and integrating. The first term on the left-hand side 
becomes dn/dt , and the right-hand side becomes the right-hand side of Eq. (2). 
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To obtain the remaining term we note that the local three-momentum squared 
is |p| 2 = gijp'pi (t.e., the physical components of the momenta are p x = Rip 1 , 
etc.), so the second term on the right-hand side is 


~g 

(2tt)3 


/ H [H lP l + H 2P : l + H 3 P l] d -£ 

= ~ 3 ( 2 *) 3(^1 + ^+^)/ 


d 3 p 

nr 


= (H 1 + H 2 + Hz)ti — 3Hn , 


(52) 


where we used the isotropy of / in the first two steps and integrated by parts in 
the third step. In doing so, we recover Eq. (2). Thus, although the form of the 
Liouville operator in the Bianchi I model differs from that in the FRW model, the 
Boltzmann equation for the evolution of the number density of a particle species 
is the same. 

The crucial assumption made above is that the particle species is in kinetic 
equilibrium. Earlier than the time of freeze out annihilations are occurring 
rapidly (r^NN > iT), and they serve to maintain both kinetic and chemical 
equilibrium. In addition, if the species (X) annihilates into relativistic particles, 
then by crossing symmetry XX's can elastically scatter with particles in the 
thermal bath with a similar cross section. The relativistic particles in the ther- 
mal bath are always more abundant than X’s— especially when 1 > 1 so that 
Y E q 1 — and so these elastic scattering processes will serve to keep XX ’s in 

kinetic equilibrium even after chemical equilibrium ceases to be maintained (i.e., 
after freeze out). 
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FIGURE CAPTIONS 


1. Plot of the Hubble parameter in a shear-dominated Universe (H) and in the 
standard model (H s t( j) as a function of temperature T. Also plotted is the 
annihilation rate r^NN a P^ticle species of mass m. The annihilation 
rate becomes equal to the expansion rate at a temperature Tf. 

2. Plot of the abundance Y of a Majorana neutrino of mass m = 60 GeV as 
a function of x = m/T in a radiation- dominated Universe (x 6 = 0), and in 
shear- dominated models with x e = 1300 and x e = 13000. 
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